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The canonical formalism in classical theory of QCD is constructed on a space-like hypersurface. 
The Poisson bracket on the space-like hypersurface is defined and it plays an important role to 
describe every algebraic relation in the canonical formalism into Lorentz covariant form. Surface 
integrals are introduced as alternatives of field equations for quarks, gluons, and Faddeev-Popov 
ghosts. It is shown that deformations of the space-like hypersurface for surface integrals are gen- 
erated by the interaction term of QCD Hamiltonian density. By converting the Poisson bracket 

Oon the space-like hypersurface to four-dimensional commutator, we can pass over to QCD in the 
Heisenberg picture without spoiling the explicit Lorentz covariance. 

(N ; 

£S| . One of the problems in the ordinary canonical method of quantization in quantum chromodynamics (QCD) is in the 
^ ■ classical procedure blinding the Lorentz covariance. The definition of the canonical conjugate momentum depends 
on the frame of reference. The Hamiltonian is linked to a definite choice of the time. Further, the commutation 
relation of two field variables is defined at the same time. From the relativistic point of view, the canonical formalism 
on the basis of a particular frame of reference spoils the explicit Lorentz covariance in the classical description of 
QCD far from field quantization. For this reason the canonical method of quantization in QCD have been regarded 
to be inferior to path-integral method. In the path-integral method Lorentz covariance is manifest from the outset. 
This is an impressive method to construct the S-matrix in QCD with covariant manner. Although the canonical 
~*^3 method is superior to the path-integral method in presenting clearly a particle interpretation to quarks, gluons, and 
Faddeev-Popov ghosts, the lack of explicit Lorentz covariance turns our interest to the path-integral method. This 
paper presents a canonical formalism in the classical theory of QCD with coherently Lorentz covariant form. The main 
" subject is to construct the canonical formalism on a space-like hypersurface. It is shown that the canonical method of 
quatization in QCD should be applied on the space-like hypersurface by converting classical fields to quantized field 
operators in the Heisenberg picture. The generalization of our formulation will be discussed in separate papers. 
The total QCD Lagrangian density with £ = 1 gauge is taken as 



+ Jflp [W*(x) {g» v d p W ip {x) - d v W^{x))] + (0"£V)) Dfc\ x ), (1) 



with 



G* (x) = fyWfrx) - d v Wh x) + gf ijk Wi( x )W*(x), 



Df = S ab d^ + gf acb W^ x ), 



where ip is the SU C (3) triplet quark, are gluon fields, c a and c a are Faddeev-Popov ghosts and anti-ghosts, A 1 
are the Gell-Mann matrices, and / yfe are the structure constants of SU C (3). With the Lagrangian density (El), the 
principle of the least action gives field equations 
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<9 M [5^iy Ip (x) + g/^W^^W^)] - gf mk W k (x)G npv (x) 
+gf ab [d p c a (x)}c b (x) - ^(,) 7 p A>(,) = 0, 



ab b . 

D abp duc\x) = 0. 



d»Dfc\ x ) = 0, 



(2) 

(3) 
(4) 
(5) 



In QCD all field variables appeared in the Lagrangian density are field operators. However, we classicize the gluon 
field operators into Grassmann even wave functions, and the quark field and Faddeev-Popov ghost operators into 
Grassmann odd wave functions to look into only the classical aspect of QCD. 

The canonical formalism needs the canonical variables. We will adopt every classical field and its conjugate mo- 
mentum on a space-like hypersurface as the canonical variables. Let us choose conjugate momenta on the space-like 
hypersurface a as 



n^-(x)=n^(x) 



d R C\. 



d [d^(x)} 
d R C[x\ 

d\d»wiv)] 



= -n p (x)d,W lv (x) - gr> k (n p { x)W>(x))W k »(x), 



d R C\x\ 

Ii a c (x)=n^x) - = n"(x)9 M c°(x), 



m(s)=n„(s) 



9[9 M c a (x)] 

d R C[x\ 
d[d^(x)] 



= -n^(x)^c°(x) - gf acb (n»(x)W°(x)) c 6 (x), 



(6) 
(7) 
(8) 
(9) 



where n M (x) is a unit normal (n M (a;)n AI (x) = 1) at x p to the surface passing through x p . The superscript " R " means 
that we have adopted the right-differentiation convention. If we make the space-like hypersurface to be global flat 
surface, the conjugate momenta @-@ reduce to the ordinary conjugate momenta. Thus our conjugate momenta are 
natural extension of the conventional ones defined at a particular time. 

The total QCD Hamiltonian density is given by the Legendre transformation on a: 



H[x\ = Il 4 ,(x)d n ^(x) + TL%(x)d n Wi(x) + Ii a c (x)d n c a (x) + Ii a £ (x)d n c a » - £[x] 



with 



H [x] = - II^(x) (n • 7) 7 jU 5^(x) - iU^(x) (n ■ 7) j v ip(x) 

- ^n w ^(x)n w ^(x) + 1 (d t ^w w (x)) (d?wi{x) 
+ n^)H>) - (d t ^ a (x)) (d?c a (x)) , 



(10) 



where 



igli^(x)(n ■ ^X'^Whx) 



f lm (n ■ W l )W Wv{x )W mv {x) - \r k Whx)W k (x) {d»W iv {x) - d?W^(x)) 



I 9 jijk jilm 



(n ■ W k )(n ■ W l )Wl(x)W mi "(x) + ^Wfa)Wj! \x)W lM (x)W m * \x) 



gf acb K(x) (n ■ W c ) c\x) - gf acb ($„<?*(*)) W c »(x)c b (x), 



d n = n ll {x)d IJ ', d tfl = {g^u - n^n^x))^. 

The differential operators d n and dtu are the directional derivatives in the direction of normal and tangent to er, 
respectively. With the help of the total canonical energy- momentum tensor T^x], we can express (fill) in the form 



H[x] = n p (x)n v (x)Tu, v [x\. 
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Thus we can identify TH,o[x] as the kinetic term of the Hamiltonian density, and Wintb] as the interaction term of the 
Hamiltonian density. The surface integral of TL on a should be regarded as the total QCD Hamiltonian: 

H = I d^(x)n v (x)T^[x] = j dSH, 

J a J (J 

where dS^ix) = n fl (x)dT,. We get Hamilton's equations by taking the variation of H on a: 



d R n 
d R n , 



d R n 
d R n 



+ d n i> = 0, 



d R n „ d R n 

dtp, 



d n U iX 



= 0, 



dip 

d R n 



QTTi d ^ 

UU IVA 



d [dtiii>] 

d R n 



d R n d R n 



d R n „ o R n 



d n Ilj, = 0, 



d n W lX = 0, 



du a c 
d R n 



d R n 



dc a 



d t ^ 7 —— + d n n a c = o, 



d n c a = o, 



d R n d R n 

~d^~ ttJ, d[d tl ,c a ] 



d n n a e = o. 



(11) 

(12) 
(13) 
(14) 



They give field equations (p) (p|) (excepting a half set of equations for gluons leading to a trivial identity = 0), so 



Hamilton's equations (|Tl )— ( 14 
Hamilton's equations (lj)- 



are equivalent to the Euler-Lagrange equations. 
14|) and the Klein-Gordon equation associated with the invariant delta function ((□- 



m )A{x) = ; OD(cc) = 0) can be combined to give the following surface integrals 

(15) 
(16) 
(17) 
(18) 
(19) 

Note that the world-point x is not necessarily on the surface c, and space-time coordinates on a are represented by 
x" . We find that the integrals (|l5|)-(^9|) satisfy free field equations: 

(17 • d — m) ip(x; <j) = 0, 
□ W iX (x;a) = 0, 
□ c a (x;a) = 0, 
Uc a {x-a) =0. 

Thus the surface integrals (|15V-(|l9|) are equivalent to asymptotic fields. 

The bubble differentiation [l] -||( a differentiation with respect to a ) of the integrals (^5|)-(|l^) around a world-point 



ip(x 


*) = 


- f d£(i 7 • d + m)A(x - x")j I^(x"), 


4>(x 


a) = 


- / dSn^{x")(-i-f ■ d + m)A(x - x"), 

J a 


W iX (x 


*) = 


[ dS[d[[D{x - x") ■ W lX (x") + D(x - x")U^(x")], 

J a 


c a (x 


<?) = 


f dS[d[[D{x - x") ■ c a {x") + D(x - x")Tl a E {x")}, 

J a 


c a (x 


*) = 


f dS[d[[D{x - x") ■ c a (x") - D(x - x")W c {x")}. 



reads 



8ij)(x; a) 
8<j{x) 
5W lX (x;a) 
da{x') 

Sc a (x; a) 
Sa(x') 

Sc a (x;a) 
Sa(x') 



= -g(i<y ■ d + m)A(x - x'^X^x'Wfa'), 

= -d'^Dix - x') ■ gf jk Wl(x')W kx (x') +D(x - x')- gf jk G jX " (x>)Wv(*') 

+D{x - x') ■ gf lcb [d' x c c (x')]c b (x') - D{x - x') ■ ig^x'^X^x'), 
= -gf acb d^D(x - x') ■ W>')cV), 

= gf acb D(x - x')WUx')d'^(x'). 



(20) 

(21) 
(22) 

(23) 
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The right-hand side of equations (p0|)-(23) includes only the interaction term, so the integrals ([D^-(|l9|) are not 
conserved because of the interaction between fields. However, equations (E0T)-(23) are invariant under the canonical 
transformations. The invariance is easily verified to represent equations (p0[)— ( 25 ) in terms of Hamilton's equations 
and the Klein-Gordon equation of the invariant delta function. 

If we impose the Lorentz condition d' p ' W lp (x") = on a as a boundary condition, the supplementary condition of 
the W ip (x;a) becomes 



d p W tp ( X] a) 



dE p ( x ")[d"^D(x - x") ■ gP k W^x")W kp (x") 

-D(x-x")-gfv k G^(,a>")Wv{*") 
- D(x - x") ■ gf cb [d' x c c {x")}c b (x") 
+ D{x - x") ■ ig^(x")-y p X l ilj(x")} w 0. 



(24) 



The weak equality symbol is to emphasize that the left-hand side is numerically restricted to be zero but does not 
identically vanish in the phase space. 

Now, we define a Poisson bracket on the space-like hypersurface a. Any two Lorentz covariant quantities F[a] and 
G[a] being the function of canonical field variables and momenta on a have a Poisson bracket which we shall denote 
by [F[a],G[a]] c , defined by 



Ja \6(j)A(z) Sn A (z) Sn A {z) 6<l>A(z) 



(25) 



The \A\ is the number of factor associated with the classical field 4>A- take for Grassmann even function, and 1 for 
Grassmann odd function. If F[a\ is given by 



F[a] = / dE"(aO fu(4>A(x),d^t4>A(x),IL A (x),d IJ ,tn.A(x)) 

J a 

= / dT, Ti^A.d^A^A.d^A) 
J a 



with dY, v f v — dYi{n- /) = dT,J-, the functional derivative on a reads 

SF[a] d R T „ d R T 



b A (z) d<f) A (z) tfJ, ~d[d tl± (f>A{z)Y 



5F[<j] d R T n d R T 



sn A (z) dn A (z) ^ d[d tl j,iiA(z)} ' 

Our Poisson bracket (|25| ) is a generalization of the ordinary equal-time Poisson bracket to be Lorentz covariant. All 
algebraic relations of the ordinary Poisson bracket hold for the Poisson bracket on a. Here we will call the Poisson 
bracket on a the four-dimensional Poisson bracket. 

Substituting (|l5|)-(|l8|) into (p5[), we get the four-dimensional Poisson bracket relations of tp(x; a), (x; a), c a (x; a), 
and c a (x; a): 

[ip a (x;a),tp l3 (y;a)] c = (ij ■ d + m) Q/3 A(x - y), (26) 
[W;(x;a),Wi(y;a)] c = -6* g„ v D(x - y), (27) 
[c a (x;a),<?(y;a)] c = 5 ab D(x-y). (28) 

If there are no interaction, ip(x; a), W*{x; a), c a (x; a), and c a (x; a) are conserved and do not depend on a. It follows 
that 

[ip a (x),iJ l3 (y)] c = (i~fd + m) al3 A(x-y), 
[Wi(x),W>(y)] c =-6» 9lMU D(x-y), 
[c a (x),c b (y)] c =5 ab D(x-y). 
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These four-dimensional Poisson bracket relations correspond to the four-dimensional commutation relations for free 
fields in QCD. 

Equations (E0h-(E3) are also expressible in terms of the four-dimensional Poisson bracket: 



Sip(x; a) 
5a(x') 
SW iX {x; a) 
da(x') 
Sc a (x;a) 

Sa(x') 
Sc a {x;a) 
Sa(x') 



- [ip(x;a),nint[x']] c , 
[W iX (x;a),H int [x']] c , 

- [c a (x;a),n in t[x']] c , 

- [c a {x;a),n int lx']] c . 



(29) 
(30) 
(31) 
(32) 



These equations show that the displacement of ip(x; a), W^(x; a), c a (x; a) and c a (x; a) under the deformation of a is 
equal to that of the surface integrals under the canonical transformation generated by the interaction Hamiltonian 
density 7i; n t at a world-point x' on a. 

To pass over to QCD in the Heisenberg picture, we shall make all canonical variables into operators and convert all 
the four-dimensional Poisson bracket relations to four-dimensional commutation relations: 

(four-dimensional Poisson bracket) — > (iU) 1 (four-dimensional commutator). 

In this way we can pass over to QCD without spoiling the Lorentz covariance. A remarkable feature of our prescription 
in quantization is that equations (29)-(|32]) are converted to the Yang-Feldman equations [Q. According to the 
procedure given by Yang and Feldman, the S-matrix of the Heisenberg picture in QCD must be identified as 



S = ■■ 



1-i I WintM d 4 x' 

CT-l 



/■CI ^ 

1-i / Wi„t[*'] dV 

J an 



where 7ii n t ty\ is the operator corresponding to the interaction Hamiltonian density 7ij n t [x' 
The supplementary condition (|24|) will be converted to 



d p W tp (x;a)- / cK p{ x")[d"^D{x ~ x") ■ gr k Wf 1 (x")W kp (x") 



D(x-x")-gf^ k G^(x")W k (x") 
D(x-x")-gf icb [d'P<?(x")]c\x») 



+ D(x - x") ■ igtp(x")'j p X l ip(x")\ (+) |* >= 0, 
where the symbol (+) means the positive-frequency part, and I'l' > is a state vector in the Heisenberg picture. 
I would like to thank Professor M. J. Hayashi for careful reading of manuscript. 
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